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Abstract. We present XMM-Newton observations of four BL Lac objects (1H1219+301, H1426+428, Markarian 501 and 
PKS 0548-322), which have been found with past X-ray missions to contain evidence of broad soft X-ray absorption features. 
Observations with the high resolution Reflection Grating Spectrometers on XMM-Newton provide the best chance yet of inves- 
tigating these features. No broad absorption features are observed in any of the objects. Neither do we find convincing evidence 
for narrow emission and absorption lines in the RGS spectra. We discuss the history of observations of broad absorption features 
in these four objects, finding that if the features exist then they must be transient, and that - given the frequency of reported 
observations of them - we can rule out the existence of transient broad absorption features in these objects at 93% confidence. 

Key words. Galaxies: active - BL Lacertae objects: individual (1H1219+301, H1426+428, Markarian 501, PKS 0548-322) - 
X-rays: galaxies — techniques: spectroscopic 



1. Introduction 

BL Lacs are well-known for their featureless spectra, all the 
way from radio to gamma rays. These smooth continua are 
thought to be generated by synchrotron self-compton emission 
from jets of material streaming out of the active nucleus to- 
wards us. Canizares & Kruper ( 1984 1, however, discovered a 
soft X-ray absorption feature in the Einstein Objective Grating 
Spectrometer (OGS) spectrum of the BL Lac PKS 2155-304, 
and since that time, similar features have been detected in sev- 
eral other objects: 1H1219+301, Ma rkarian 501, PKS 0548- 
322, Markarian 421 (Madejski et al. ITWfl and H1426+428 
(Madejski et al. 1992 1. This absorption, interpreted in each 
case as being due to ionised gas intrinsic to the BL Lac, is per- 
haps the only insight we have into the immediate environment 
of these enigmatic sources. The dispersive spectrometers on- 
board XMM-Newton and Chandra, with their combination of 
high resolution and large effective area, provide the best op- 
portunity yet to investigate and understand these features. To 
this end, XMM-Newton observed the BL Lacs 1H1219+301, 
H1426+428, Markarian 501 and PKS 0548-322 with the aim of 
obtaining high resolution soft X-ray spectra with the Reflection 
Grating Spectrometers (RGS; den Herder et al 2001 1, and high 
signal-to-noise CCD resolution spectra from EPIC (Striider et 
al. |2001). This paper presents the X-ray spectra of these ob- 
jects, and discusses the evidence for spectral traces of the BL 
Lac environment. Table^gives their coordinates, cosmological 
redshifts and Galactic neutral absorbing columns. 

Various authors have searched for soft X-ray ionised ab- 
sorption in the four BL Lacs discussed in this paper. Madejski 



et al. ( 1991 ) found evidence for an absorption notch (trough) 
in the Einstein SSS spectrum of 1H1219+301; this feature ap- 
peared at 0.56 +0 02 keV (all energies in this Section are quoted 
in the observed frame unless otherwise stated) and had a width 
of about 100 eV. It has apparently not been seen since with 
other missions. 

A soft X-ray absorption feature in the BBXRT spectrum of 
H1426+428 was reported by Madejski et al. i 19921 1. and anal- 
ysed in greater detail by Sambruna et al. ( 1997). These latter 
authors also find indications of soft X-ray ionised absorption 
in spectra from ROSAT and ASCA. The BBXRT feature was 
best described by a broad gaussian absorption line at ~0.66 keV 
with a width less than or equal to 43000 km s -1 . The ROSAT 
and ASCA data were fitted together as they had similar flux lev- 
els. The fitting produced indications of two absorption edges at 
0.15 keV and 0.52 keV. 

The first evidence for an absorption feature in 
Markarian 501 came from the Einstein SSS spectrum; 
Urry et al. (1986 1 found an absorption trough at ~0.57 keV 
with a width of 80 eV (these were the same as the parameters 
obtained for PKS 2155-304). The data were reanalysed by 
Madejski et al. I199H . who modelled the absoption with 
a notch at 0.56 + 0.01 keV, of width 240 eV and covering 
factor 0.52. A later study by Fink et al. fi99~fl of ROSAT data 
revealed the possibility of an absorption line at 0.37 keV 
(the width was fixed at 100 eV), although the statistical 
significance was too low to prove the existence of the line. 

Madejski et al. ( 1991 ) detected an absorption trough at 0.57 
± 0.02 keV in the Einstein SSS spectrum of PKS 0548-322. 
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More recently, Sambruna et al. (1998 1 found an absorption 
edge at ~0.66 keV in ASCA data; this feature could alterna- 
tively be explained by a notch at 0.82 keV with width 100 eV 
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Table 1. The coordinates, redshifts and Galactic neutral Hydrogen columns of the sources. 



Target 


Right Ascension 


Declination 


Redshift a 


Galactic column b 


1H1219+301 


12h 21m 21.941s 


+30d 10m 37.11s 


0.182 


1.78° 


H1426+428 


14h 28m 32.6s 


+42d 40m 21s 


0.129 


1.36 d 


Markarian 501 


16h 53m 52.2s 


+39d 45m 37s 


0.03366 


1.73° 


PKS 0548-322 


5h 50m 40.771s 


-32d 16m 17.76s 


0.069 


2.51" 



a redshift as listed by the NASA Extragalactic Database (NED) 

b Galactic absorbing column in 10 20 cirT 2 

c Elvis etal. H989l 

d Costamante et al. 120011 



(fixed) and covering factor 0.12. They re-analysed the Einstein 
SSS data, obtaining an edge at 0.487 keV, and found evidence 
of two edges in ROSAT data. One of the ROSAT edges (at 
around 1 keV) was probably instrumental, and the other, at 
0.553 keV, may either have been intrinsic or due to residual 
calibration problems. 

The usual interpretation advanced for these features is that 
they originate in highly ionised oxygen in the environment 
of the BL Lac. The earliest detected feature amongst the ob- 
jects described here, the absorption trough in Markarian 501, 
was simply described as being due to highly ionised oxygen 
(Urry et al. [T986l . Madejski et al. fTWil . following the lead of 
Canizares & Kruper ( 1984 1, interpreted their notches as being 
due to O vin Lyo- resonance absorption. Sambruna et al. ( 1997 1 
found their -0.66 keV feature to be consistent with absorption 
edges of O v-O vn, and marginally consistent with O vm Ly<x 
In their ROSAT and ASCA fits, no physical meaning could be 
determined for the 0.15 keV edge, whereas that at 0.52 keV 
was identified with the K-edge of Om. Fink et al. (1991 1 give 
no physical interpretation for their 0.37 keV line. The absorp- 
tion feature seen by Sambruna et al. (1998) in PKS 0548-322 
has two possible explanations: if it is an edge, it corresponds to 
O vi, and if it is a notch it would correspond to Fe xix-Fe xx or 
Ne vn. Sambruna et al. (I1997l and ll998i also fit the absorption 
with photoionised warm absorber models. The appearance of 
features at different energies, when the sources are at different 
flux levels, has also been used as a diagnostic. In H1426+428, 
the features were found to be at lower energies when the source 
flux was lower, consistent with the predictions of warm ab- 
sorber models (Sambruna et al. 1997 1. The opposite has been 
found to be true (lower energy features at higher flux levels), 
however, in PKS 0548-322 (Sambruna et al. ll998l 

Seyfert galaxies have the best understood warm absorbers 
of all the classes of AGN. Since the advent of high resolution 
X-ray spectrometers, it has been known that narrow absorp- 
tion and emission lines, rather than deep absorption edges, are 
the most common characteristic of ionised gas in the environ- 
ment of this variety of AGN (e.g. Kaastra et al. 2000 1. It would 
be interesting to see if this is also true of BL Lacs, and par- 
ticularly whether the low resolution of previous spectrometers 
has been hiding a wealth of information in these narrow lines. 
Since the launch of Chandra, weak intrinsic narrow absorp- 
tion lines of Ovn Hea have been found in Chandra LETGS 
spectra of Markarian 421 (Nicastro et al. 2000 2003 1, and nar- 



row absorption lines due to the local Warm-Hot Intergalactic 
Medium have also been found in PKS 2155-304 (Nicastro et 
al. l2002> and Markarian 421 (Nicastro et al. l2003> as observed 
by Chandra, and also in XMM-Newton RGS spectra (Paerels 
et al. 2003). There have been no reports so far of the broad ab- 
sorption features found in the past, even though PKS 2155-304 
and Markarian 421 have been observed many times as calibra- 
tion sources by both XMM-Newton and Chandra. 

2. Observations and data analysis 

The dates, observing modes and exposure times of the observa- 
tions are given in Table|2] The EPIC-pn data were processed us- 
ing the XMM-Newton SAS (Science Analysis System) Version 
5.4, using the epproc task, and spectra and lightcurves were 
extracted from a source region of 40", using both single and 
double events, with background subtraction being performed 
using regions of the same size. Time intervals affected by pro- 
ton flaring were filtered out of the event lists used for spec- 
tral extraction; the medium filter exposures of Markarian 501 
were disregarded due to a high proton background throughout 
in both of them. Response matrices and ancillary response files 
were generated using the SAS tasks rmfgen and arfgen re- 
spectively. It was found, using epatplot, that the minimum 
useful energy for spectral fitting is ~0.4 keV. There are no pn 
data for PKS 0548-322. 

The EPIC-MOS data were extracted with emproc under 
SAS V5.4. The imaging mode MOS data for all objects were 
found to be piled-up, and so were not used further in the analy- 
sis. Timing Uncompressed (TU) mode data were available for 
H1426+428 and Markarian 501; these did not suffer from pile- 
up, although the response matrices were found to have an in- 
correct energy scale in the region of the oxygen K-edge at ~ 
0.5 keV which biased the spectral fitting. Because of this, and 
the existence of better-quality pn data, the TU MOS data were 
not used in the final spectral fits. 

RGS1 and RGS2 data were extracted using rgsproc under 
SAS V5.4, which carries out background subtraction using re- 
gions spatially adjacent to the source. Time intervals affected 
by proton flaring were filtered out. For each source, RGS 1 and 
RGS2 data (both first and second order spectra) were combined 
channel by channel to achieve the best possible signal to noise. 
The RGS spectra were binned into groups of three channels, 
approximately equivalent to the spectral resolution. Both the 
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RGS and pn X-ray spectra were analysed using xspec, with 
the search for detailed spectral features being carried out in 
SPEX 2.00 (Kaastra et al. l2002at . Since PKS 0548-322 has 
no usable EPIC data, the RGS was used to obtain an X-ray 
lightcurve with evselect, performing background subtraction 
using the average background lightcurve obtained from regions 
above and below the source in the cross-dispersion direction. 

3. Lightcurves 

The 0.4-10 keV pn lightcurves of 1H1219+301, H1426+428 
and Markarian 501 are constant to within ~ 5%, and the RGS 
lightcurve of PKS 0548-322 is on average constant to within 
~ 10% (the signal-to-noise is poorer for the RGS lightcurve); 
these lightcurves are plotted in Fig^ Our spectral fitting is thus 
unlikely to be much affected by any spectral variability in the 
sources. 

4. EPIC-pn spectra 

The 0.4-10 keV pn spectra were fitted with power-law models 
with fixed neutral absorption due to our Galaxy. The ratios of 
the data to these models are shown in Fig.|2] Some spectral cur- 
vature is present in each case, but there is no evidence for deep 
warm absorption - or excess neutral absorption - in any of the 
sources; any fit residuals in the soft band are all within the 10% 
calibration uncertainties which exist below 2 keV. Motivated 
by the observed spectral curvature, we fitted broken power-law 
models (again with fixed Galactic neutral absorption) to the 
spectra. The parameters of both broken power-law and simple 
power-law fits are listed in Table [3] and the broken power-law 
fits are shown superimposed on the data in Fig.|2] 

5. RGS spectra 

None of the RGS spectra show evidence of broad absorption 
features. The parameters of power-law and broken power-law 
continuum fits, with neutral absorption due to the ISM of our 
Galaxy, are given in Table [3] The statistical significance of ap- 
parent narrow emission and absorption features in the spectra 
was assessed by fitting a gaussian against a local continuum (as 
done for NGC 7469 in Blustin et al. l2003l . The results of this, 
alongside the respective RGS spectra, are plotted in Figs. [3J- 
[6] If the features were purely due to statistical noise, then in 
an RGS spectrum of about 1000 data points, one would expect 
to see about three features significant at 3<x. Any 4<x features 
would be likely to be real. 

Figs. [5] and |5] show that the apparent features in 
1H1219+301 and Markarian 501 are consistent with statisti- 
cal noise. H1426+428 (Fig.0) and PKS 0548-322 (Fig.® each 
have one narrow absorption feature significant at 4<x or greater. 
The H1426+428 feature has an equivalent width of 40 + 10 mA 
at 1 1 .55 A (observed frame), which is 10.23 A in the source rest 
frame, and could be identified as Ne x LyyS. The corresponding 
Lyo- is not present though, making the identification unlikely. 
The observed line could conceivably be Ne x Lye blueshifted 
by ~ 47000 km s -1 , if it was associated with the BL Lac jet, 
but no further lines are found which could give support to this 



pattern. The 4<x absorption line in PKS 0548-322 appears at 
21.82 A in the observed frame (20.41 A in the rest frame), and 
has an equivalent width of 70 + 20 mA at 21.82 A. It has no 
obvious interpretation in the rest frame - the nearest likely tran- 
sition is Ovu Hea, at a blueshift of ~ 17000 km s _1 . Again, 
there are no further features to reinforce this interpretation. 

H1426+428 and PKS 0548-322, respectively, have five and 
ten features significant at three sigma. In each case, only one 
of the features is in emission. This is not what is expected from 
a random distribution, and suggests that the features may be 
real. They are not due to bad channels in the spectrum, and the 
vast majority do not seem to be related to known systematic 
errors (e.g. anomalies in the response matrix at CCD gaps). We 
have, however, been unable to find any consistent or plausible 
identification or interpretation for any of the features. 

6. Discussion and conclusions 

Our high resolution and CCD resolution spectra showed no 
sign of deep broad ionised absorption features, as had been 
claimed to be present in spectra from previous missions. Two 
RGS spectra, those of PKS 0548-322 and H1426+428, have 
weak indications of narrow line absorption, although we are 
not able to find convincing identifications for any of these ap- 
parent features, and so cannot confirm that they are real. 

So where have the deep intrinsic absorption features gone? 
Fig-Qcontains plots of the observed 2-10 keV X-ray flux ver- 
sus time, over the past 25 years, of the four sources discussed 
in this paper. Those observations where features have been re- 
ported are marked. Only published observations by missions 
that had sufficient soft X-ray effective area and spectral res- 
olution to detect the features are included (Einstein, ROSAT, 
BBXRT, ASCA, BeppoSAX and XMM-Newton; there have 
been no Chandra publications on these objects). It has been 
claimed (e.g. Madejski et al. I199U that intrinsic absorption 
could be detected in EXOSAT spectra, by means of a soft X- 
ray spectral flattening, even though this mission had no energy 
resolution in the band where the features occur; we find no re- 
ports in the literature that EXOSAT has ever actually detected 
such absorption for any of our four objects, so EXOSAT ob- 
servations are not included. ROSAT's energy resolution in the 
relevant energy band is extremely low, but since ROSAT spec- 
tra of H1426+428 and Markarian 501 have shown possible ev- 
idence of ionised absorption (Sambruna et al. 119971 Sambruna 
et al. 1998 Fink et al. 1991}, we include flux measurements 
from this mission. The fluxes, when not quoted directly in the 
2-10 keV band in the original papers, were calculated from 
the original fit parameters using either PIMMS or XSPEC, and 
the derived values are given in Table 0] alongside the observa- 
tion dates and references. The 2-10 keV band was chosen as 
it is the one quoted most frequently in the literature (thus min- 
imising the chance of introducing errors in the extrapolation to 
a new band), and also the range least affected by the possible 
confusion between power-law slope and excess neutral absorp- 
tion in the soft band. 

The data in Fig.^suggest that the features - if they exist - 
are only visible some of the time. However, absorption could 
have been present but not reported (or not detectable) during 
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Table 2. Observation dates, X-ray instrument modes and exposure times for the BL Lac observations reported here. 



Object 


observation date 


EPIC-MOS1 3 


EPIC-MOS2 b 


EPIC-pn c 


RGSl d 


RGS2 C 


1H1219+301 


2001-06-11 


29304 (SW M) 


29304 (FF M) 


28523 (SW M) 


29854 


29854 


HI 426+428 


2001-06-16 


4044 (TU M) 
60407 (LW M) 


4044 (TU M) 
60357 (FF M) 


66239 (SW M) 


67442 


67454 



Markarian 501 2002-07-12 
2002-07-14 



6802 (TUT) 8699 (FF M) 6910 (SWT) 10232 10227 
1012 (TU M) 1720 (SW M) 

7177 (TUM) 11877 (FFM) 6910 (SWT) 12654 12647 
3857 (TUT) 4201 (SW M) 



PKS 0548-322 2001-10-03 



47303 (LW M) 47277 (FF M) (SW M) 



48152 48161 



a Exposure time in seconds for MOS1. The operating mode and filter are given in brackets (SW = Small Window, TU = Timing 

Uncompressed, LW = Large Window; M = Medium filter, T = Thick filter). 
b Exposure time in seconds for MOS2. The operating mode and filter are given in brackets (FF = Full Frame, TU = Timing Uncompressed, 

LW = Large Window; M = Medium filter). 
Exposure time in seconds for pn. The operating mode and filter are given in brackets (SW = Small Window; M = Medium filter, T = Thick 

filter). 

d Exposure time in seconds for RGS 1 ; all exposures were in Spectroscopy mode. 
c Exposure time in seconds for RGS2; all exposures were in Spectroscopy mode. 
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Fig. 1. The 0.4-10 keV pn lightcurves of 1H1219+301, H1426+428 and Markarian 501 (from orbits 474 and 475), all in 300 s 
time bins, and the RGS lightcurve of PKS 0548-322 (0.35-2.5 keV; 600 s time bins). 
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Table 3. Fit parameters for powerlaw and neutral absorption, and broken power-law with neutral absorption, fits to the EPIC-pn 
(0.4-10 keV) and RGS (0.35-2.5 keV) spectra of the BL Lacs objects observed. The 2-10 keV flux is given for the preferred 
model for each object. 



Object 


Instrument 


N H a 






r 2 d 


N e 


y 2 f 

A red 


F2-10 8 




1 TTI ^1 A . T A 1 

1H1219+301 


pn 


1.78 


2.485 ± 0.006 






2.23 ± 0.08 


1.38 (714) 








pn 


1.78 h 


2.40 ± 0.02 


1.4 ± 0.2 


2.61 ± 0.02 


2.29 ± 0.02 


0.99 (712) 


2.64 ± 


0.02 




RGS 


1.78 h 


r\ f\r\ r\ f\r-\ 

2.09 ± 0.02 






2.00 ± 0.02 


1.08 (1038) 








KLro 


I. 10 


l./O ± U.Uo 


n ~iz _i_ n n^ 
U. / J ± U.U3 


Z.44 ± U.U / 


Z.J ± U.l 


U.oy (lUio) 






H1426+428 


pn 


1.36 1 


1.870 ± 0.003 






0.899 ± 0.002 


1.61 (607) 








pn 


1.36 


1 mo i a nnn 

1.922 ± 0.009 


1.5 ± 0.2 


1 1 C i A f it \( \ 

1.816 ± 0.009 


r\ ofll i A AA1 

0.891 ± 0.003 


1.24 (605) 


2.90 ± 


0.01 




RGS 


1.36' 


1.68 ±0.02 






0.778 ± 0.006 


1.19 (996) 








RGS 


1.36' 


1.47 ±0.07 


0.79 ± 0.08 


1.88 ±0.05 


0.89 ± 0.05 


1.05 (994) 






Markarian 501 


pn 


1.73 h 


2.308 ± 0.005 






2.639 ± 0.008 


1.17 (653) 








pn 


1.73 h 


2.34 ± 0.01 


1.9 ±0.4 


2.23 ± 0.03 


2.62 ± 0.01 


1.03 (651) 


4.48 ± 


0.02 




RGS 


1.73 h 


2.00 ± 0.02 






2.27 ± 0.02 


1.01 (1038) 








RGS 


1.73 h 


1.70 ±0.06 


0.71 ± 0.04 


2.26 ± 0.05 


2.8 ±0.1 


0.84 (1036) 






PKS 0548-322 


RGS 


2.51 1 


1.74 ± 0.02 






0.942 ± 0.008 


1.28 (1038) 








RGS 


2.51' 


1.38 ±0.05 


0.74 ± 0.02 


2.01 ± 0.04 


1.19 ±0.04 


1.03 (1036) 


2.47 ± 


0.08 



a neutral absorbing column in 10 20 cirr 2 

b slope of soft power-law 

c break energy in keV 

d slope of hard power-law 

c power-law normalisation in 10~ 2 photons keV" 1 cm" 2 s _1 at 1 keV 

f reduced x~ of the total fit; the degrees of freedom are given in brackets 

g 2-10 keV flux of preferred model in 10~ 11 erg cirr 2 s~' 

h fixed at Galactic Hydrogen column (Elvis et al. 1989 1 

1 fixed at Galactic Hydrogen column (Costamante et al. l200H 



some of these observations, especially the earlier ones. Because 
the XMM-Newton spectra are of much higher quality than any 
of the previous data, and no broad features are observed in any 
of the XMM-Newton spectra, we can state categorically that if 
the features observed with past missions were real, they must 
be a transient phenomenon. 

In Fig. no clear patterns emerge in terms of feature visi- 
bility at different flux levels. In H1426+428, a feature of some 
sort was observed for at least three years, but is no longer seen 
eight years after it was first detected. In Markarian 501, the 
absorption appears in a space of a few months and has disap- 
peared by the time of the next observation four years hence. We 
can use the detection rate of the broad absorption features in 
previous observations to determine the probability of observ- 
ing no features with XMM-Newton. The fraction of previous 
observations which did not yield a detection of a broad ab- 
sorption feature is 75%, 25%, 66% and 57% for 1H1219+301, 
H1426+428, Markarian 501 and PKS 0548-322 respectively. 
With these detection rates, the probability of observing no 
broad absorption features in any of these objects with XMM- 
Newton is only 7%. 

If we take these results at face value, then we can rule out 
the existence of transient deep, broad intrinsic absorption fea- 
tures in these sources at 93% confidence. If a further compa- 
rable set of observations fail to detect any features, then their 
existence will be ruled out at 99.5% confidence. Time - and 
future observations with the high quality instruments on board 



XMM-Newton and Chandra - will lay the question of broad 
absorption features in BL Lacs to rest. 

Low spectral resolution, low signal-to-noise and calibra- 
tion uncertainties must have contributed to previous conclu- 
sions about these features. The major source of the apparent de- 
tection of intrinsic soft X-ray absorption features was, though, 
probably the use of a single power-law to fit the X-ray spec- 
tral continuum of BL Lac objects. These continua were fitted 
at higher energies and extrapolated over the soft band, and the 
resulting deficit of photons in the soft band was interpreted as 
evidence for absorption. As our observations show, there is sig- 
nificant spectral curvature in our four sources. This is also seen 
in other recent observations; for example, Fossati et al. (p000l 
have shown that the X-ray continuum of Markarian 42 1 is best 
described by a gradually steepening power-law with the photon 
index increasing towards higher energies. Such spectral curva- 
ture, when fitted with a simple power-law, would give rise to 
soft X-ray residuals which could plausibly be interpreted as ab- 
sorption. 

So, BL Lac objects do not seem likely to yield much knowl- 
edge of the immediate environment of their central engines 
through X-ray absorption features. A more fruitful search for 
such evidence will probably involve other classes of radio 
galaxy in which we are not looking directly down the axis of 
the jet; Broad Line Radio Galaxies with Seyfert-like nuclei are 
an ideal candidate for such studies. 



6 



A. J. Blustin et al.: Intrinsic absorbers in BL Lac objects: the XMM-Newton view 




Biergy (ksV) ^Tgy 

Fig. 2. Left panels: ratios of the pn spectra of 1H1219+301, H1426+428 and Markarian 501 to fitted power-laws with fixed 
Galactic neutral absorption (plotted in the observed frame). Right panels: the pn spectra of 1H1219+301, H1426+428 and 
Markarian 501 with superimposed best-fit broken power-law models (with fixed Galactic neutral absorption; plotted in the ob- 
served frame). In all plots, the spectra are grouped to a minimum of 40, 350 and 100 counts per bin respectively. 
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Fig. 3. The statistical significance of narrow absorption and emission features in the RGS spectra of 1H1219+301 (top), and the 
RGS spectrum itself with best-fit broken power-law model (red) superimposed (bottom), plotted in the source rest frame. 
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Fig. 4. The statistical significance of narrow absorption and emission features in the RGS spectra of H1426+428 (top), and the 
RGS spectrum itself with best-fit broken power-law model (red) superimposed (bottom), plotted in the source rest frame. 
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Fig. 5. The statistical significance of narrow absorption and emission features in the RGS spectra of Markarian 501 (top), and the 
RGS spectrum itself with best-fit broken power-law model (red) superimposed (bottom), plotted in the source rest frame. 
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Fig. 6. The statistical significance of narrow absorption and emission features in the RGS spectra of PKS 0548-322 (top), and the 
RGS spectrum itself with best-fit broken power-law model (red) superimposed (bottom), plotted in the source rest frame. 
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Fig. 7. 2-10 keV fluxes of 1H1219+301, H1426+428, Markarian 501 and PKS 0548-322 observed with Einstein, ROSAT, 
BBXRT, ASCA, BeppoSAX and XMM-Newton. Observations where intrinsic absorption was reported are circled. 
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Table 4. The history of the 2-10 keV fluxes of 1H1219+301, H1426+428, Markarian 501 and PKS 0548-322 as observed by 
past X-ray missions. 



Object 


Mission 


Date 


F(2-io) a 


Intrinsic ab- 


Reference for flux value 










sorption 




1H1219+301 


XMM-Newton 


11/06/01 


2.6 


no 


Tabled 




BeppoSAX 


13/07/99 


1.5 


no 


Costamante et al. 2001 




ROSAT 


18/06/92 


2.7 


no 


Lamer et al.119961 




Einstein 


05/06/79 


1.1 


yes b 


Urry et all 19861 




Einstein 


07/12/78 


2.5 


no 


Urry et all 19861 


H1426+428 


XMM-Newton 


16/06/01 


2.9 


no 


Table[3] 




BeppoSAX 


08/02/99 


2.0 


no 


Costamante et al.12001 1 




ASCA 


06/02/94 


1.6 


yes 


Sambruna et al. 1997 




ROSAT 


24/01/92 


2.3 


yes 


Sambruna et al. 1997 




BBXRT 


09/12/90 


3.5 


yes 


Sambruna et al.119971 


Markarian 501 


XMM 


13/07/02 


4.5 


no 


Tabled 




ASCA 


26/03/96 


6.6 


no 


Kubo et al. [19981 




ROSAT 


13/08/91 


2.2 


yes 


Fink et al.119911 




ROSAT 


25/02/91 


1.3 


no 


Lamer et al. 11 9961 




Einstein 


22/08/79 


2.7 


yes 


Ciliegi et al.119931 




Einstein 


01/03/79 


5.8 


no 


Ciliegi et al.119931 




Einstein 


25/01/79 


1.8 


no 


Ciliegi et al.|1993l 


PKS 0548-322 


XMM 


03/10/01 


2.5 


no 


Table|3l 




BeppoSAX 


07/04/99 


1.5 


no 


Costamante et al.120011 




BeppoSAX 


26/02/99 


1.8 


no 


Costamante et al.120011 




BeppoSAX 


20/02/99 


2.3 


no 


Costamante et al.120011 




ASCA 


30/10/93 


2.8 


yes 


Sambruna & Mushotzky|1998l 




ROSAT 


06/03/92 


3.8 


yes c 


Lamer et al.ll996l 




Einstein 


06/04/79 


4.0 


yes 


Sambruna & MushotzkvfT998l 




Einstein 


10/03/79 


4.5 


no 


Ciliegi et alPEM! 



a 2-10 keV flux in 10" 11 erg enr 2 s" 1 

b Madeiski et al. l!99"Tl 

c Sambruna & Mushotzky |1998l 



